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This work explores a fundamental dynamical structure for a wide range of many-body quantum 
systems under periodic driving. Generically, in the thermodynamic limit, such systems are known 
to heat up to infinite temperature states in the long-time limit irrespective of dynamical details, 
which kills all the specific properties of the system. In the present study, instead of considering 
infinitely long-time scale, we aim to provide a general framework to understand the long but finite 
time behavior, namely the transient dynamics. In our analysis, we focus on the Floquet-Magnus 
(FM) expansion that gives a formal expression of the effective Hamiltonian on the system. Although 
in general the full series expansion is not convergent in the thermodynamics limit, we give a clear 
relationship between the FM expansion and the transient dynamics. More precisely, we rigorously 
show that a truncated version of the FM expansion accurately describes the exact dynamics for a 
certain time-scale. Our theory reveals an experimental time-scale for which non-trivial dynamical 
phenomena can be reliably observed. We discuss several dynamical phenomena, such as the effect of 
small integrability breaking, efficient numerical simulation of periodically driven systems, dynamical 
localization and thermalization. Especially on thermalization, we discuss a generic scenario on the 
prethermalization phenomenon in periodically driven systems. 


I. INTRODUCTION 
A. Physical background 

Periodically driven quantum dynamics has recently 
attracted much attention in experimental as well as 
theoretical studies [1-5] as it offers a promising way 
for exploring novel quantum phenomena which would 
be difficult or impossible to observe otherwise. Al¬ 
though the instantaneous Hamiltonian at each time 
step is very simple, dynamical behavior can be highly 
nontrivial. Remarkable dynamical phenomena include 
dynamical localization [6-10], coherent destruction of 
tunneling [4, 11, 12], localization-delocalization transi¬ 
tion [13-17], and dynamical phase transitions [18-22]. 
Moreover, recent experimental development has rapidly 
opened new possibilities to control quantum systems 
under periodic driving, e.g., in the context of quantum 
transport [5, 23-25], quantum topological phases [26- 
32] and detections of the Majorana Fermion [33, 34] 
and the Higgs mode in condensed matter [35]. 

One of main subjects in driven quantum many-body 
systems is to understand the thermodynamical struc¬ 
ture of steady states. With a few exceptions [10, 18- 
22], recent studies mostly focus on driving simple (of¬ 
ten non-interacting particles) Hamiltonians. However, 
the integrability-breaking terms unavoidably exist in 
the realistic experimental conditions. When one looks 
at long-time behavior, even small integrability-breaking 
terms are relevant to the dynamics and cause signifi¬ 
cant effects on the final steady states. This provides 
the motivation to understand the long-time behavior 
of ‘generic’ many-body systems under periodic driv¬ 


ing. Long-time behaviors in driven non-integrable sys¬ 
tems are in general very complicated and cannot be 
captured with simple techniques such as the rotating- 
wave approximation and the transfer matrix technique 
of Landau-Zener transitions [6] etc., and hence one is 
obliged to rely on numerical calculations. 

Recently, true steady states in the long-time limit 
have been intensively studied using large scale numer¬ 
ical calculations [10, 16, 36]. In the long-time limit, 
periodically driven many-body systems are in general 
expected to heat up to infinite temperature. This is a 
consequence from the analogy of eigenstate thermaliza¬ 
tion hypothesis (ETH) in non-driven many-body sys¬ 
tems [37]. The ETH implies that each energy eigen¬ 
state of Hamiltonian is indistinguishable from the mi- 
crocanonical ensemble with the same energy. In period¬ 
ically driven systems, the energy is no longer conserved 
and hence the extension of ETH to the driven case in¬ 
dicates that the steady state is a state of infinite tem¬ 
perature (i.e., the completely random state). Thus, due 
to the heating effect, any information reflected from the 
system is invisible in the infinite-time scale. 

However, the experiments on many-body quantum 
systems do not focus on the long-time limit; rather, they 
are interested in the transient dynamics for the experi¬ 
mental time scale. In this time scale, the heating pro¬ 
cess cannot necessarily occur, and hence a critical task 
that follows is to clarify the time-scale during which one 
can observe transient behavior that can show nontrivial 
phenomena. So far, most studies on transient dynami¬ 
cal properties are based on numerical calculations with 
phenomenological arguments. 

In this paper, we establish a new framework that de- 
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FIG. 1. In the thermodynamic limit, the FM-expansion (1) 
diverges in general (See Fig. 2 for the numerical demonstra¬ 
tion). Nevertheless, as we show, the truncated expression 
(2) at some order no accurately describes the finite-time be¬ 
havior of the periodic dynamics. 


scribes transient dynamics for a wide range of many- 
body systems. Our analysis is based on the Floquet- 
Magnus theory. We show that as long as we consider a 
finite-time scale, the time-evolution of the system is ap¬ 
proximately governed by a well-defined effective Hamil¬ 
tonian. We will give a prescription to obtain the ef¬ 
fective Hamiltonian and reveal the time scale for which 
the system is governed by this Hamiltonian. With the 
technique, we also discuss several dynamical phenom¬ 
ena relevant to the transient dynamics including ther- 
malization process in driven systems. 


B. Outline of the basic framework 

In our analysis, the Floquet theory plays a fundamen¬ 
tal role [38] in defining the effective Hamiltonian that 
describes stroboscopic time-evolution for every period. 
This effective Hamiltonian is referred to as the Floquet 
Hamiltonian H-p and is defined as 

g-iffpT _ 

with T the time-ordering operator, where the Hamilto¬ 
nian is periodic in time: H{t) = H{t + T). Through¬ 
out the paper, we set H = 1 and to = 0) and we 
do not consider the degree of freedom on the Floquet 
gauge [1, 31, 39, 40]. The Floquet Hamiltonian has 
full information on the thermodynamic properties of 
the driven systems, and hence finding this is one of the 
central problems in periodically driven systems. How¬ 
ever, the time-ordering operation in the unitary time- 
evolution is, in general, difficult to analyze, and hence 
finding the Floquet Hamiltonian is a highly nontrivial 
problem [1, 39-44]. 

The Floquet-Magnus (FM) expansion [45, 46] is 
known to give a formal expression of the Floquet Hamil¬ 
tonian as follows 

OO 

Hp = T^fln (FM expansion), (1) 

n—0 

where explicit forms of the terms {i4n})][Lo 9'^'® given 
by Eq. (6) in Subsection HB. In this paper, we also 


consider the truncated expansion up to the n-th order, 
which is defined as 

n 

iTp"^ ;= T™Q.m (Truncated FM expansion). (2) 

m—0 

The FM expansion (1) is useful especially for the high- 
frequency limit in finite-size systems, where the higher- 
order contribution is negligible. On the other hand, 
in the case of finite frequencies, the higher-order con¬ 
tribution should always be considered. However, it 
has been recognized that, in general, the FM expan¬ 
sion is not convergent series expansion [1, 46]. The 
divergence problem in the FM expansion has been a 
long-standing obstacle in analyzing periodically driven 
systems, and hence our understanding of the Floquet 
Hamiltonian for many-body systems is so far extremely 
limited [10, 16, 36]. 

The divergence of the FM expansion is expected to 
have an interesting physical meaning [10, 13, 15-17, 36]. 
In periodically driven systems, even the energy is not a 
conserved quantity. Hence, the ergodic time-evolution 
implies that the entire phase space can be covered 
by the time evolution from an arbitrary initial state. 
References [13, 36] have numerically shown the signa¬ 
ture of ergodic dynamics in the time-evolution opera¬ 
tor showing that it is described by a random matrix. 
Based on these observations, one anticipates a deep 
connection between the divergence of the FM expan¬ 
sion and the quantum ergodicity [10, 13, 15-17, 36]. 
This connection indicates that the divergence is usually 
unavoidable because generic non-integrable many-body 
systems are believed to be ergodic in the thermody¬ 
namic limit [36, 37, 47-51]. 

Even though the FM expansion is divergent in gen¬ 
eral, we empirically expect that the finite truncation of 
the FM expansion (2) can give useful information on the 
transient dynamics. We will show that the FM expan¬ 
sion has an optimum order to approximate the transient 
dynamics even if it is divergent in the limit of n ^ oo; 
the figure 1 schematically shows the typical behavior of 
the FM expansion. The validity of the finite trunca¬ 
tion of the FM expansion has been studied in several 
specific cases such as the Friedrichs model [52] and the 
NMR of solids [53]. Our purpose in this paper is to 
generalize these results and to give a rigorous relation¬ 
ship between the FM expansion and general properties 
of transient quantum dynamics. 

In the present study, we consider a generic isolated 
quantum system on a lattice, where we do not impose 
any assumption such as the integrability or symmetry 
conditions. We only assume that the interaction cou¬ 
plings between spins (particles) are few-body in Hamil¬ 
tonian. This class of Hamiltonians has, thus far, been 
employed to study on quantum information oriented 
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problems [54-58]. Here, we apply the techniques de¬ 
veloped for such problems to the Floquet theory, and 
rigorously show the condition under which the trun¬ 
cated expression (2) can be used and estimate its accu¬ 
racy, towards the main aim of this study. As shown 
in Table I, we will show that by choosing the trun¬ 
cated Hamiltonian H^°'^ with ng = 0{uj), the error 

of II”*^|| decays exponentially as w in¬ 
creases. This fact leads to the notion of quasi-stationary 
states for long but finite time-scale; the system first re¬ 
laxes to an intermediate state which is characterized by 
the effective Hamiltonian H^°\ and then reaches the 
final steady state (i.e., the infinite temperature state) 
very slowly. We will see that such a quasi-stationary 
state is maintained up to an exponentially long time 
scale of 

The present paper is organized as follows. In Sec. H, 
we show the setup of systems and our main results on 
the generic properties of the transient dynamics. As 
shown in Sec. HI, our results are related to various facts 
in fundamental physics. In Sec. IV, we summarize our 
results and mention areas for further development. Fi¬ 
nally, Section V is devoted to proving the main theo¬ 
rems. 


II. SETUP AND MAIN RESULTS 


where hx is an operator acting on a spin subset X. 
More explicitly, it can be given in the form of 


H — ^ ^ hiSi -f ^ ] Jil ^i 2 A ^ ^ 'Ai ,22 ,^3 ^ 

i ii,i2 * 1 , 22,*3 


^i2 


^ik ! 


where {si} are operator bases on the Ah spin; for exam¬ 
ple, it can be given by the Pauli matrices for (l/2)-spin 
systems, namely {s^} = {trf, erf, erf}. This definition 
of the Hamiltonian encompasses almost all interesting 
quantum many-body systems with short-range interac¬ 
tions such as the the XY model [59], the Heisenberg 
model [60, 61] and the AKLT model [62], as well as 
models with long-range interactions such as the Lipkin- 
Meshcov-Glick model [63]. Typically, we have k = 2 
(i.e., two-body interaction), but several exceptions ex¬ 
ist such as the cluster-Ising model [64] (fc = 3), the toric 
code model [65] (fc = 4) and the string-net model [66] 
{k = 6). Although we do not treat fermionic systems 
explicitly, our discussions can be also applied to local 
fermionic systems because they can be mapped into lo¬ 
cal spin systems [67, 68]. 

We then introduce a parameter J as a local interac¬ 
tion strength (or one-particle energy) of the system: 


In this section, we first introduce the setup of our sys¬ 
tem and show the convergence of the Floquet-Magnus 
expansion up to a finite order with several numerical 
demonstrations. We then give our main theorems on 
the transient dynamics. 

A. Generic few-body Hamiltonians 

We consider a spin system of finite volume with each 
spin having a d-dimensional Hilbert space, and we label 
each spin by i = 1,2,... N. We denote the set of all 
spins by A = {1,2,... N}, a partial set of sites by X, 
and the cardinality of X, that is, the number of sites 
contained in X, by |A| (e.g. X = {H, * 2 , ■ • ■, *|x|})- We 
also denote the complementary subsets of X by X^; that 
is, X 0 X^ = A. Note that the spins are not necessarily 
located next to each other on the lattice. 

We consider a general periodically driven system on 
an arbitrary lattice, the Hamiltonian of which is given 
by H(t) = Ho-\-V(t), where Hq is the static Hamiltonian 
without the driving potential and V{t) = V{t-\- T) de¬ 
notes the driving Hamiltonian with the period T (with 
the frequency uj := l/T). Throughout the paper, we 
consider the system which is governed by a generic few- 
body Hamiltonian, that is, the Hamiltonian contains at 
most fc-body interactions with finite k: 

Ho = X! X! (^) 

|x|<fc |.^|<fc 


X! Vie A, (4) 

X:X3i 

where || • • • || is the operator norm and J^x xbi denotes 
the summation with respect to the supports containing 
the spin i. Note that we thereby obtain J2x x 3 i \\^x\\ < 
J and condition indicates 

that the energy change due to one spin is bounded by 
a finite value J. Throughout the paper, we use the 
notations 

Vo ■■= Y. f r A:=2fcJ, 

|X|<fc ^ 

0{uj) = const • Y and 0{T) = const • XT, (5) 
A 

where VqT denotes the driving amplitude in one pe¬ 
riod. In our theory, the parameter A characterizes typ¬ 
ical properties of the system. 


B. Floquet-Magnus (FM) expansion 

A possible method to calculate the Floquet Hamilto¬ 
nian Hp is to expand it with respect to the period T as 
in Eq. (1), namely Hp = A"fin- The FM expan¬ 

sion gives each of the terms {H„}fTo as follows [69]: 


4 


^77, - 


1 

{n+iy 

1 


jn^n+1 


^ n\ 

(7 

[ dtn+l...[ dt2 [ dti[H{t^(^ri+l)),[H{tcr(n)),- ■ ■ ,[H{t„:^2)),H{t^(^i^)] . . .]] , 

Jo Jo Jo 


( 6 ) 


where a is the permutation and 0(a) := X]r=i 0(a(i + 1) — a(i)), with O(-) the usual step function. For example, 
the first three terms in the expansion read (See [70] for higher-order terms) 


ffo(T)=^^ i7(ti)dti, ffi(r) = ^^ dhJ^'dt2[Hih),Hit2)], 

^ 2 ( 7 ) = —^^^ J dti J dt2 j dt3([77(ti), [iJ(t2),77(t3)]]-I-[i7(t3), [i7(t2),i7(ti)]]). 


Without loss of generality, we can put jJ'v(t)dt = 0 
and, thus, flo = Hq. The FM expansion is considered 
a useful tool to treat a periodically driven system when 
the period T of the driving is sufficiently small. 

From a general discussion, we can evaluate the upper 
bound of ||T”n„|| as 

where ^ is a universal constant [46, 71-74]. Then, the 
convergence of the FM expansion is ensured only for 
the case of /J' l|i7(t)l|dt < and it is not satisfied for a 
macroscopic system (x N), unless the period T 

scales with the total system size N. Thus, this simple 
estimation indicates the divergence in the thermody¬ 
namic limit. 

When the Hamiltonian is given by a few-body opera¬ 
tor (3) with the condition (4), much stronger bound for 
exists: 

Lemma 1. For each of the terms with n > 1, 

oy \n _ 

||f^n||< (7) 

(n + 1}^ 

where A and Vq are defined in Eq. (5). The proof is 
given in Appendix C. 

The bound of Eq. (7) implies that the FM expan¬ 
sion is convergent up to n ~ (AT)“^ in the sense of 
l|n„llT” > for n < However, for 

n > (Ar)“^, the FM expansion begins to diverge. In 
the next subsection, we show that even if the FM ex¬ 
pansion might diverge for n —> 00 , the truncation (2) 
characterizes the quantum dynamics with great accu¬ 
racy. 


C. Numerical demonstration in anisotropic 
Heisenberg chain 

Before showing our main theorems, we would explain 
our results in a visual way. For the purpose, in Fig. 2, we 
give numerical calculations of the FM expansion. Here, 



iog(ii4"'ii) 


r = 0.5 T = f).4 


(b) Logarithmic plot of ||t?p”^| 



FIG. 2. The logarithmic plots of (a) norms of each term 
in FM expansion, (b) norms of the truncated Flo- 

quet Hamiltonians (c) norms of ^ 

with respect to the order n. We consider the Heisenberg 
model with a time-periodic field for three cases of the pe¬ 
riod T = 0.2,0.3,0.4,0.5, which correspond to the purple, 
blue, orange, and red curves, respectively. 


we consider an 8-spin anisotropic Heisenberg chain with 
a time-periodic field: 




for 0 < t < T with H(t) = H(t + T), where {a^}^=x,y,z 
are the Pauli matrices and we assume the periodic 
boundary condition, namely ug = tJi. We calculate the 
FM expansion for four cases of T = 0.2,0.3, 0.4,0.5 up 
to the order n = 40. 

We plot the norms of the expansion terms 
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(Fig. 2 (a)) and the truncated Floquet Hamiltonians 
(Fig. 2 (b)). The calculation implies that the 
FM expansion is neither conditionally convergent nor 
absolutely convergent. Thus, for T > 0.3, the Mag¬ 
nus expansion completely breaks down for n ^ 1, 
while the Magnus expansion converges up to a small 
n. This is consistent to the Lemma 1. Relevant to 
them, we demonstrate the accuracy of the truncated 
FM expansion in Fig. 2 (c), where we plot the error 
Ijg-iffpr _ g-i^Tp T’li respect to the truncated or¬ 
der n. We can see that the qualitative behavior is syn¬ 
chronized with that of ||T"n„|| in Fig. 2 (a). In the 
following, we show the analytical reasonings to these 
observations. 


D. Transient dynamics for generic few-body 
Hamiltonians 

We shall prove the following theorem on the valid¬ 
ity of the truncated Floquet Hamiltonian for the time 
evolution: 

Theorem 1. Consider Hq and V(t) = V(t + T) be 
few-body operators with the local interaction strength 
J, respectively, and assume that T satisfies T < 1/(4A). 
Then, the time evolution under the Hamiltonian H(t) = 
Hq + V(t) is close to that under the truncated Floquet 
Hamiltonian H^°'^ = with 


in the sense that 

||g-iffpT _ ^ 6Vbr2-"“, (8) 

where [-J is the floor function. See Eq. (5) for the defi¬ 
nition of A and Vq. We give the proof in Section V C. 

Note that the right hand side of (8) is exponentially 
small for w, as VoTe~^^‘^\ Although the above theorem 
compares the time evolutions only for one period, it 
immediately follows that for t = mT with m a positive 
integer: 

||g-*ffpmT _ < QVomT2-^°. 

Thus, this theorem indicates that the exact time evolu¬ 
tion is well approximated by the time evolution under 
the ng-th order truncated Floquet Hamiltonian up to 
time t « Even if the Floquet-Magnus expan¬ 

sion is not convergent and the system might eventually 
heat up to infinite temperature, the transient dynam¬ 
ics of the system is governed by the truncated Floquet 
Hamiltonian for an exponentially long time with respect 
to the frequency of external driving. 

In practice, we are often interested in the dynamical 
behavior by with n less than rig; indeed, we mostly 
treat the case of n = 0, namely the average Hamilto¬ 
nian. By using Eq. (7), we can also know how close 
the exact time evolution is to that under the n-th order 
truncated Hamiltonian: 


Corollary 1. Under the conditions in Theorem 1, we 
have 

||g-iffpT _ e-*4"’T|| < 6VbT2-"« -h (9) 

for an arbitrary n < uq, where is defined in Eq. (7). 
We show the proof in Appendix A. 

This corollary rigorously shows that the Floquet- 
Magnus expansion is, at least, an asymptotic expansion. 

E. Transient dynamics for short-range few-body 
Hamiltonians 

In Theorem 1, we have shown that the total uni¬ 
tary operator is approximately described by 

° However, the inequality (8) is meaningful 
only when VgT2“”“ < 1, or uj > logVg; because Vg is 
proportional to the number of spins subjected to driv¬ 
ing, the condition is not satisfied for finite frequencies 
under global driving (Vg oc N) in the thermodynamic 
limit of fV —>■ oo. In the following, we improve the con¬ 
dition by imposing a stronger restriction to the Hamil¬ 
tonian, namely the assumption of the short-range inter¬ 
action. 

We here focus on the time evolution of a local region 
L instead of the total system A; let p be an arbitrary ini¬ 
tial state and consider its reduced density matrix in the 
region L, namely pL = tric(p), where tric(- • •) denotes 
the partial trace operation with respect to the spins in 
L^. We now define 

PL{mT) = tric(e-*"^^pe*^^^), 

(mT) = tr^c (10) 

and obtain the upper bound of \\pL{mT) — p^^°\mT)\\i 
with II • 111 the trace norm. Under the condition of the 
short-range interaction, the error between pL{mT) and 
p^^°\mT) is exponentially small for w up to t « 
as long as |T| < 

In considering short-range interacting systems, we 
have to define the structure of the system explicitly 
(e.g., the square lattice) [75]. We now define a set of 
the bonds Af,, i.e. pairs of spins {ziU 2 }, {* 2 , * 3 }, {* 2 , * 5 } 
and so on. The form of A^ decides the structure of the 
lattice. Based on this definition, we define the distance 
dist(A, Y) as the shortest-path length which one needs 
to connect the two partial sets X and Y. 

In this subsection, we introduce the following addi¬ 
tional assumption to the Hamiltonian: 

\\hx\\<F{r) for Vi e A, 

X:X9i,diam(X)>r 

where diam(A) := sup^j j}ex dist(z, j) and the function 
F{r) determines how the interactions decays as the spa¬ 
tial distance increases. In this case, we can prove the 
Lieb-Robinson bound for arbitrary operators Ox and 
Oy [75-79]: 

\\[Ox{t),Oy]\\ < g{l,t)min{\X\,\Y\)\\Ox\\ ■ ||Oy||, 

( 11 ) 
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where I = dist(X, Y) and Q{1, t) is characterized by the 
form of the interaction decay F{r). Note that is 

a monotonically increasing function with respect to t. 
The Lieb-Robinson bound characterizes the non-locality 
of operator due to the time-evolution; that is, we can 
know how fast an operator Ox spreads away from the 
region X after a short time. 

Using the above notations, we can obtain the follow¬ 
ing theorem for short-range interacting Hamiltonians: 

Theorem 2. Let us consider a Z?-dimensional system 
which satisfies (11). Then, the exact time evolution of 
PL{mT) by is close to the approximate time 

evolution p^^°\mT) by °^mT 

<12 J|L|mT2-”«/2 -h 2\L\mg{lo, mT) (12) 

with Iq = const where H-Hi denotes 

the trace norm. We give the proof in Section VD. 

For example, in the case where the interaction decays 
exponentially {F(r) ^ e“’’), we know [75-78] 

(13) 



FIG. 3. Schematic picture of Theorem 3 for H^\ Let us 
consider an initial state \'4><e) which is in a superposition of 
energies below E (blue curve). Then, after a time of t = mT, 
the energy distribution still localizes in the range [0, E + 
VoO{T’^'^^) + beyond which the distribution 

decays exponentially (red curve). 

the eigenspaces of (n < no) with energies larger 

than E. Then, for an arbitrary initial state \iP<e) which 
is in a superposition of energies below E (see Fig. 3), 
we define the probability to absorb energy AE after 
m periods as ||n^^^^p,|V’<£;(t))|P, where t = mT and 
\'4’<Eit)) = Now, Theorem 1 is immedi¬ 

ately followed by 

l|ng+ABlV^<ijW)f < i6tVo2--of . (14) 


with c, V and ^ some constants depending on A and 
E{r). Then, we have g{lQ,mT) = cexp[—(6®^“^/^^ — 
umT)/^]. Hence, by applying the Lieb-Robinson 
bound (13) to the inequality (12), the second term is 
negligibly small in comparison with the first term as 
long as mT < Thus, in this time-scale, the 

error between pL{mT) and p^^°\mT) is as small as 

On the other hand, when the interaction decays poly- 
nomially {E{r) ~ l/r“), the Lieb-Robinson bound is 
given by [75, 79] 

g{l,t) = c— for D < a < 2D, 

g{l, t) = -b c for a>2D 

with {c, c',u,^, 7 } constants depending on A and F{r). 
Thus, from the inequality (12), in the former case, the 
error \\pL{mT) — p^^°\mT)\\i is exponentially small 
g-O(w) time t ft! O{io), while in the lat¬ 

ter case, the error is exponentially small up to the time 
t ft This way, if the interaction decays rapidly 

enough (a > 2D), we can ensure that the truncated FM 
expansion gives a good approximation as long as we look 
at the local region of the system (|L| < We em¬ 

phasize that Theorem 2 does not depend on the driving 
amplitude Vq and is meaningful even for global driving 
(Vo oc N) in the thermodynamic limit of —>■ oo. 


for n = uq. In this estimation, the probability does not 
depend on the amplitude of AE; however, we expect 
that the probability of energy absorption should depend 
on the energy amplitude AE. Indeed, we can prove the 
following much stronger statement than (14): 

Theorem 3. When we assume T < r := 1/(8A) with 
A := 6k^J, the following inequality holds 


IIH 


>E+AE\'^<E{t)) II 


< exp 


2t (^-AE + UoT”+^VV„+i -b 22tXVo2-^°/^'j 


with t = mT, where Wn is defined by 1U„ := 
2(4A/3)"'n!/(n-bl)^. We show the proof in Appendix E. 

From Theorem 3, we can identify the following two 
points (Fig. 3): 

1. After a time t = mT, the system can absorb en¬ 
ergy at most AE ft VoO{T^~^^)+tXVoe~‘^^‘^'>; that 
is, the energy absorption rate is exponentially sup¬ 
pressed for the frequency. 

2. The probability for large energy absorption decays 
exponentially or even faster as AE increases. 

On the point 1., Ref. [80] has recently obtained the simi¬ 
lar results in the linear response regime, which has been 
generalized in our recent paper [81]. 


F. Probability of energy absorption 

We here discuss an energy absorption by the peri¬ 
odic driving. Let H^^ be a projection operator onto 


III. RELATIONS TO SEVERAL PHENOMENA 

From our results, we can identify several properties 
which universally appear in periodically driven systems. 
We here relate them to the existing results in literatures. 












7 


1. Comparison with experimental time scale. 

Experimentally, the control parameter of the system is 
often given hy hw/J (See references in [1]). Remember 
that the parameter J characterizes the one-spin energy. 
From Theorem 1, the accuracy of the approximation is 
ensured as long as 

^ < 2V(i6AT) ^ (^5) 

with k = 2 (i.e., two-body interaction), where we ex¬ 
plicitly describe the Planck constant h. Hence, the re¬ 
liable time of the approximation sensitively depends on 
the coefficient of hio/J, which is now given by log 2/64. 
When we consider an optical field as the periodic driving 
(w « Hz), our present estimation shows that, 

as long as fujj/J ~ 10^-10^, the lifetime mT (or m/w) 
is comparable to the experimental time scale. This con¬ 
dition still leaves room for improvement as we presently 
consider the most general setup in estimating the coef¬ 
ficient in (15). 

2. Thermalization. Based on our results, we give a 
general scenario for the thermalization process in peri¬ 
odically driven systems. For arbitrary n with n < no, 
as shown in the inequality (9), the time evolution of 
the n-th order truncated Floquet Hamiltonian approx¬ 
imates the exact time evolution up to t « 0(w"+^). 
For larger t, the time evolution is not approximated by 

However, up to t « H^°'^ still gives a good 

approximation for the time evolution. As long as H^°'^ 
is ergodic, long-time unitary time evolution under 
would result in an equilibrium state [37, 48-50], that is, 

{tp{t)\0\^{t)) ~ tr(Op[;(/>)), (16) 

where O is an arbitrary local operator and pmc ^ is the 
density matrix corresponding to the microcanonical en¬ 
semble of Because the difference between 

and i7p"°^ is tiny for small T, the microcanonical ensem¬ 
ble of i7p”^ would approximately equal to that of H^°\ 

Pmc — Pmc^ for any n < uq. We therefore expect that 
for sufficiently small T, we have 

(V’(i)|O|'0(t)) ~ tr(Op[;)/) for Vn<no (17) 

as long as t < Although |||V'(i)) — '*!'*/')II 

is not small for t ^ the system will reach a 

quasi-stationary state represented by the microcanon- 
ical ensemble of Hp and this quasi-stationary state 
persists at least up to t Thus, even the 0-th 

order truncated Floquet Hamiltonian would well 
describe a quasi-stationary state of a periodically driven 
system up to an exponentially long time with respect to 
w. 

This characterizes a kind of prethermalization in the 
sense that the system first relaxes to a quasi-thermal 
state Pmc \ which is far from the final equilibrium state. 


namely the infinite temperature state. The mechanism 
of this prethermalization originated from the approxi¬ 
mate freezing of the energy absorption and is qualita¬ 
tively different from the established ones [82-86]. 

3. Dissipative systems with periodic driving. 

Our theorem is also applicable to an open system in 
which a small driven system is in contact with a large 
heat bath, as long as the total system is governed by 
a few-body Hamiltonian and the state vector of the to¬ 
tal system obeys the Schrodinger equation. According 
to the discussion in thermalization, the total system 
including the heat bath goes to the microcanonical en¬ 
semble, which implies that the system of interest reaches 
a quasi-stationary state described by the canonical en¬ 
semble of Hp ^ (the Floquet-Gibbs state [51, 87]) with 
the inverse temperature /3 identical to that of the ther¬ 
mal bath and this quasi-stationary state persists at least 
up to t At first glance, it seems contradictory 

to the fact that the long-time asymptotic state of the 
system of interest is in general not given by the Floquet- 
Gibbs state even for large w [87-91] . To resolve it, we 
should remember that the system of interest stays in 
the Floquet-Gibbs state for a fixed-time scale, given by 

gOM^ 

Driven open systems are often treated by the Floquet- 
Born-Markov master equation [51, 87], which is valid in 
the van Hove limit [92] (or the weak-coupling limit); as 
shown below, we cannot use the Floquet-Born-Markov 
master equation to grasp the finite-time behavior of 
the system. If the interaction Hamiltonian between the 
small system and the thermal bath is denoted by XHj, 
the limit of A —)■ 0 and t —> oo with X^t held fixed is 
called the van Hove limit. Therefore, in the van Hove 
limit for a fixed value of w, the time scale of t « 
is always infinitesimally shorter in comparison with 
the relaxation time, and thus the Floquet-Born-Markov 
master equation cannot capture the quasi-stability of 
the Floquet-Gibbs state, which has the lifetime of 
This necessitates us to treat the dissipative system be¬ 
yond the van-Hove limit in order to discuss the realistic 
relaxation process. This necessitates us to treat the 
dissipative system beyond the van-Hove limit in order 
to discuss the realistic relaxation process (see also the 
recent paper [93]). 

4. Numerical simulation in low-energy regime. 

Our analysis includes an implication on efficient numer¬ 
ical calculation to simulate periodically driven systems. 
Our theorem implies the existence of approximate local 
conserved quantities, which is given by Prac¬ 

tically, we often cannot obtain the high-order Magnus 
terms for large systems, whereas only a few terms of 
the series give us a lot of useful information. For exam¬ 
ple, let us choose the ground state of Hp^ as the initial 
state. Then, after a time of t < Theorem 3 en¬ 

sures that the energy is exponentially concentrated in 
the range of [Eo,Eo 0{T)Vo), where we denote the 
ground state energy of E[p^ by Eq. This fact is bene¬ 
ficial in numerically simulating the low-energy periodic 
dynamics, because the entanglement area law is ensured 



for any low-energy state in various classes of Hamilto¬ 
nians which are characterized by the heat capacity [94] ; 
thereby, we can apply the matrix product state or the 
projected entangled pair state to simulate the dynamics 
up to the time t « [95-98]. 

5. Dynamical localization. When the dynamical 

localization occurs, the energy of the system cannot ex¬ 
ceed a certain non-trivial energy bound. From Theo¬ 
rem 3, we can define as the approximately con¬ 
served energy for a long time of t « hence the dy¬ 

namical localization generally occurs for large frequen¬ 
cies, but in the approximate sense. We notice that the 
evidence of such behavior has been numerically reported 
in Ref. [10]. 

The dynamical localization often occurs in periodi¬ 
cally driven systems, whereas the mechanism is not the 
same in all of the cases; for example, the dynamical 
Anderson localization [7] occurs because a periodic kick 
effectively causes random potential in momentum space, 
the dynamical freezing or coherent destruction of tun¬ 
neling [11, 99-101], originate from the resonant inter¬ 
ference effect at particular frequencies, the many-body 
localization [13-15, 102-104] is related to the local in¬ 
dependence of spins which comes from disorder in the 
Hamiltonian, etc. 

In the present case, we can provide a qualitative ex¬ 
planation in terms of the following viewpoint. In the 
periodically driven system, the unit of energy which the 
system can absorb is roughly quantized by w, while the 
energy associated to one spin is bounded by J as in (4) . 
Thus, for the energy absorption of 0{uj), many spins 
(i.e., uj/J spins) should work synergistically in one pe¬ 
riod. However, such many-body quantum effect is ex¬ 
ponentially suppressed in a short-time evolution for the 
few-body Hamiltonians [58]. In this way, the energy 
absorption rate is bounded from above by this 

mechanism of the (quasi-) dynamical localization is dis¬ 
tinct from the existing ones. 

6. Small integrability breaking. We here con¬ 
sider the case where a small integrability breaking term 
eV is added to an integrable Hamiltonian H*{t) = 
H*{t + T), that is, H{t) = H*{t) + eV, where eV is 
time-independent. Then, the Floquet Hamiltonian Hp 
is formally expanded as 

Hf := + ei7i,F + e^H2,F + ■■■ , 

where Hp is the Floquet Hamiltonian with respect to 
which might have a non-trivial form far from 
the random matrix. For example, let us choose bilinear 
fermionic Hamiltonians (e.g., XY model) as the unper¬ 
turbed Hamiltonian H*(t). Then, the Floquet Hamilto¬ 
nian Hp is always far from the random matrix [105, 106]; 
that is, heat up to the infinite temperature never oc¬ 
curs. On the other hand, an infinitesimal perturbation 
eV can cause the heating after sufficiently long time, 
which kills all the specific properties which are origi¬ 
nated from Hp. We discuss the time-scale for which 
the integrability breaking term becomes dominant. 


In order to give the effect of eV, we further expand 
Hj^F with respect to T as 

where can be determined by the use of the FM 

expansion. Then, our result ensures that the system is 
governed by the following effective Hamiltonian 

no 

=H^ + eJ2 + en2,m + • • •) 

m—0 

as long as t < This implies the following: 

i) Up to the time t « 1/e, the system is well approx¬ 
imated by the integrable Hamiltonian Hp and de¬ 
scribed by the generalized Gibbs ensemble [105]. 

ii) The non-integrable terms + 

+ • • •) influence the dynamics in the time 
scale of 1/e ^ t jt and reduce the sys¬ 

tem to the Gibbs ensemble by H^°\ In this time 
scale, the effective Hamiltonian H^°'^ is still close 
to Hp, and hence the characteristic properties of 
the system may strongly depends on Hp. 

iii) The divergence terms of the FM expansion finally 

kill all the non-trivial structure of Hp for t ^ 
gO(i^)y'g yp |-]^g systom to the infinite 

temperature state. 


IV. SUMMARY AND FUTURE PERSPECTIVE 

In summary, we have worked on the generic tran¬ 
sient dynamical behavior of periodically driven systems. 
In the thermodynamic limit, the Floquet-Magnus ex¬ 
pansion (2) usually diverges for n —> oo, whereas it 
converges up to a certain order of 0{uj) as shown in 
Lemma 1. Our main results characterize the time-scale 
for which the truncation of the FM expansion provides 
a good approximation: ~ ° ^ with an ex¬ 

ponentially small error for the frequency. It is worth 
noting that our results can be applied to macroscopic 
systems {N —> cx)) and are highly universal in that we 
consider generic few-body Hamiltonians without spe¬ 
cific assumptions. Our results have given the first the¬ 
oretical step to understand the transient dynamics sys¬ 
tematically. 

On the other hand, our work has left several open 
problems. First, the main theorems only refer to the up¬ 
per bound on the error of approximation, and hence the 
infinite-time behavior is still an open problem. From the 
convergence condition of the FM expansion, we can en¬ 
sure that only when co > N the exact Floquet Hamilto¬ 
nian is far from the random matrix. On the other hand, 
our results indicate that for w > log TV the truncated 
FM expansion gives a very good approximation at 

least for the unitary operator, i.e., ~ 

Hence, we expect that we might obtain a looser con¬ 
dition than ui > N for the exact Floquet Hamiltonian 
to have a well-defined form. Indeed, the recent results 
support this expectation; for example, the finding by 
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D’Alessio and Rigol [36] implies that the long-time be¬ 
havior of a periodically driven lattice system qualita¬ 
tively changes at around w > 'i/N. 

Second, it would be an interesting direction to extend 
our present theory to Liouvillian dynamics in dissipa¬ 
tive driven quantum systems. The dissipation causes 
non-trivial quantum phenomena which cannot be sim¬ 
ply explained and has recently attract much attention 
not only in fundamental [107-109] but also in experi¬ 
mental aspects [110-112]. We believe that our frame¬ 
work will be applicable to generic Liouvillian dynamics 
with time periodicity; here, we can define an effective 
Liouvillian operator. The problem is not straightfor¬ 
ward because the basic assumptions of (3) and (4) are 
not trivially ensured for generic Liouvillian operators. 
Even though, several techniques applied here, such as 
the Lieb-Robinson bound, can be also applicable to var¬ 
ious classes of dissipative dynamics [113, 114]. We ex¬ 
pect that the qualitative properties are the same be¬ 
tween the effective Liouvillian operator and the Floquet 
Hamiltonian in closed systems. 

Third, related to the discussions 2 and 3 in Sec. Ill, 
it is an intriguing problem whether our scenario on 
the thermalization can be numerically (or experimen¬ 
tally) observed or not. Without a heat bath, a pe¬ 
riodically driven system is conjectured to show two- 
step relaxation process; it first goes to a prethermalized 
state which is characterized by the effective Hamiltonian 
and then relaxes to the final steady state, i.e. the 
infinite temperature state. On the other hand, in the 
presence of a heat bath, the quasi-steady state charac¬ 
terized by might have infinite lifetime, where the 

energy emission to the heat bath should be comparable 
to the energy absorption from the external driving. Al¬ 
though we will have to treat a considerable number of 
spins in the numerical demonstration, we expect that 
we can simulate the relaxation processes efficiently in 
the low-energy regime according to the discussion 4 in 
Sec. HI. 


V. METHODS 

In this section, we show the proofs of our main the¬ 
orems. We first show a decomposition of the Floquet 
operator as an important stepping stone to the proof. 
Based on it, we give outlines of the proofs and defer the 
details of the calculations to Appendix sections. 


A. Decomposition of the Floquet operator 

As a main idea for the proofs, we introduce a decom¬ 
position of the total unitary operator: 

C/p = c-^^o^UnUn-i ■ ■ ■ U 2 U 1 (18) 


rr('^o) _ —iHQTjj{'^o)jT{'i^Q) rr('^o) 7-r(^o) 

U-p — c ’ 


where we denote the unitary operators e and 

T jj^ Lp"\ respectively. In this decom¬ 
position, we first separate the driving Hamiltonian V(t) 
into N pieces, namely {Fi(t)})li, and aim to treat their 
driving effects separately. As shown below, the decom¬ 
posed unitary operators correspond to the con¬ 

tributions by {Vi{t)}^i, respectively. 

We first introduce the notations of Ag := 0, A^ := 
{1, 2,..., f} and 

Mt) ■= -VKit), 

N 

A(t):=Ho-h ^ y,(t) (20) 

j=i+l 

for i = 0,1,...,iV, where Ho{t) = H{t), HN{t) = Hq 
and Hi{t) = Vi+i{t) + Note that Vi{t) is 

the driving operator which contains the spin i but not 
contains the spins {1,2, ••• ,i — 1}. Using the above 
notations, we first decompose the total Hamiltonian 
into H{t) = Vi{t) -I- which reduces the total 

unitary operator to C/p = Uf,i(T)Ui with Up^iit) := 

7-[e-*/o andC/i = 7-[e”*-/o^. 

Here, Up^i{T) is the Floquet operator with respect to 
the Hamiltonian We can then decompose as 

H^{t) = V 2 {t) + H 2 {t) and define Uf,i{T) = Uf, 2 {T)U 2 
in the same way. By repeating this process, we de¬ 
compose the total unitary operator Up as in Eq. (18), 
where the decomposed unitary operators are given by 
UF,i-i{T) = UF^i{T)Ui such that 

U- •= T[e~*/o^ 

C/F,,(t) := r[e-*/o ■ ■ ■ U,+i (21) 

for 1 , 2 ,..., A. 

Second, we define as the truncated FM expan¬ 

sion at Tigth order with respect to the time-dependent 
Hamiltonian where Hq^°^ = and 77^“^ = 

Hq because of Ho{t) = H{t) and Hnii) = Hq, re¬ 
spectively. We then introduce the notation of := 

fj(no) _ gino) 

N 

i7f”«) = i7o+ (22) 

j=i+i 

Based on the above notations, we also decompose the 
effective unitary operator in the similar way to 

(21): 

Jj(no) _ 

U^y\t) := (23) 

for i = 1, 2,..., A. Note that the definition (23) implies 

C^f72i(7’) = 

For the convenience, we often include Hq into V{t); 
that is, we denote H{t) = Hq + V'{t) with Hq = 6 


and 
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and V'{t) = H{t). In this case, the decomposition (18) 
reads 

= (24) 

Then, the definition (20) implies that the Hamiltonian 
H'iit) — contains the spins {i + l,z + 

2 ... ,N}, and hence from (21), the unitary U- is sup¬ 
ported in the spin set {i,i + 1..., N}. 

B. Basic lemma on the decomposed unitary 
operators 

For each of the decomposed Floquet operators 
{Ui}^i, we can derive the following lemma; we utilize 
it as a central technique for the subsequent proofs. 

Lemma 2. Let and be the de¬ 

composed unitary operators given in (21) and (23). 
Then, each of the unitary operators {Ui}^^ is close 
to with no = [jg^J as follows: 

<6V,r2-"«, (25) 

where A is defined in Eq. (5) and we define as 

V.:= E ^r\\vxm<J. (26) 

X-.X3i,XnAi_i=</) ”'0 

Notice that because of Eq. (5) 

N 'Y' 

E^*=E^/ II^^WII = ^o. (27) 

i=i xeA ^ 

Proof of Lemma 2. For the proof, we start from the 
equalities C/f.*-i(T) = and = 

from which we expand Ui and 

with respect to the period T: 

OO 

= C/t ,(T)[/f,._i(T) = 

n—0 

CXD 

uino) ^ (28) 

n—0 

where the terms and are invariant for the 

scaling T —>• sT, H{f) —> H{t/s) with s a positive con¬ 
stant. 

Now, the operator U^°\t) comes from the FM ex¬ 
pansion for the Hamiltonian Hi(t). Hence, from the 
basic property of the FM expansion, when we expand 
f/F,i(T) and U^f\T) with respect to T as UF,i{T) = 

and 4J\t) = respec¬ 

tively, the expanded terms are identical up to the order 
no, namely Ui^n = for n < no- By applying the 

above equality to (28), we also obtain 

for n < no, 


which yields 

OO 

\\U.-UY^\\< Y T’”(||<i>..n|| + ||<„‘’^||). (29) 

n=no + l 

We then evaluate the upper bound of ||4)i_„|| and 
based on the explicit forms of (21) and (23) 

for Ui and uY\ respectively. After straightforward 
but rather technical calculations, we can obtain the fol¬ 
lowing inequalities for the norms of expanded terms: 

r”||4>,,„(r)|| < V,Tei/(2fc)2-n-H^ (30) 

< V*rei/(®'=)2-”+b (31) 

We show the derivations in Appendix D 1. 

By applying the inequalities (30) and (31) to (29), we 
obtain 

||[/. _ uY'^W < 2V*r2-"“(ei/^ -f e^/®) < 6V*r2-"», 

where we use fc > 1 from the definition of (3). This 
completes the proof. ■ 

C. Proof of Theorem 1 

By applying Lemma 2 to the decompositions ( 18) and 
(24), the proof is immediately followed by 

\\Uf - < \\U 1 U 2 ■■■Un- “^11 

<\\Ui- \\ + \\U2---UN-Ut°'’ ■■■U^x°'>\\ 

N N 

< Y \\Ui - uY'^W < 6T2-"« Y (32) 

i=l i=l 

By using Eq. (27) in the above inequality, we complete 
the proof. ■ 

D. Proof of Theorem 2 

For the proof, we first show the following corollary on 
the time evolution of a local operator: 

Corollary 2. Let Ol be an arbitrary operator sup¬ 
ported in a region L. Then, the exact time evolution of 
Ol{T) by Uf is close to the approximate time evolution 
by uY^-. 

\\OAT) - O^A°\t)\\ < 12J|L|r2-"«||Oi||. 

Proof of corollary 2. In order to prove the corollary, 
we consider the decomposition (24), and then we can 
obtain the same inequality as (25): 

<6JT2-"“. (33) 

We now label spins in the region L by {1,2,..., |L|}. 

Remembering that from the definition (24) 

are supported in the region we have Ol{T) = 
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Ol / OLit) 



\l Li 

OL{t) ^ OL{t,l) 


FIG. 4. Schematic picture of the Lieb-Robinson bound. We 
consider an operator Ol which is initially supported in a 
region L. after a short time, the operator Ohit) is no longer 
supported in the region L. However, we expect that the 
operator On(t) still remains in a region closed to L; that 
is, it is approximated in some region Li having distance I 
from L. The Lieb-Robinson bound gives the accuracy of 
this approximation. 


where has been defined in Eq. (11). 

We now define Sm ■= \\OLimT) — O^L°\mT)\\ and 
obtain for Sm+i 

Sm+I <\\UlOL{mT)UF - 

+ ||0L(mT)-0[”‘’^(TOT)|| 

<\\ulOL{mT,l)UF - 
+ 2\\OL{mT) - OL{mT,l)\\ + sra 
<12 J\Li\T2-'^‘> +2\L\g{I, mT) + s„„ 

where we apply Corollary 2 and (34) in the last inequal¬ 
ity. We thus obtain 

s™ < 12J|L;|rm2"”'>-b2|L|mg(Z,mT), (35) 


UlOlU'l with the notations of U'^ = ■■ - in 

the same way, we also have (T) = UOlU 

From the inequality (33), we obtain \\U'l — < 

6J|L|T2“"“ in the similar way to the inequality (32). 
We thus obtain 

\\Ol{T) - 0i”“HT)|| = Wu'^OlU'l - 
< 2||;7[ - ■ \\Ol\\ < 12J|L|r2-"“||Oi||. 

This completes the proof of Corollary 2. ■ 


Proof of Theorem 2. From the definition (10), pL^t) — 
is an Hermitian operator, and hence we can 
always find an operator Ol with ||Ol|| = 1 such that 


\\PL{t) - Pi “^(t)||i = tr Ol {pL{t) - pt°\t)) 


= tr 


p{OL{t) 


-o[”«)(t))l < ||0i(t)-0^"») 


where we use |tr(i?A)| < ||il|ji • |jA|| for arbitrary opera¬ 
tors A and B; note that ||p|ji = 1. We therefore have to 
calculate the upper bound of ||Oi(TOT) — Oi"‘’^(TOr)||. 

For the evaluation of the norm, we can apply Corol¬ 
lary 2 to ||Ol(T) — 0^l°\t)\\, whereas for \\Ol{‘2T) — 
o[”“)(2T)|| = ||l/FOL(T)t/j!: - 

the corollary cannot be applied directly; that is, the 
operator Ol{T) is no longer supported in the region L 
due to the time evolution (Fig. 4). We thus would like to 
know how accurate the operator OL(t) can be approxi¬ 
mated by an operator acting on spins which is close to 
L. For this purpose, we denote the set of spins having 
distance at most I from L hy Li, namely dist(L,i) < I 
for Vz S Li. We define the following operator 0^(1, 1) in 
the region L; which approximates the operator OL(t): 

OL{t,l) := } tric[C>i(t)] 0 lic. 

tr(lic) 


From Ref. [77] , the Lieb-Robinson bound gives the error 
of the approximation as 


(34) 


where we use Sq = 0 and g{l,mT) < 0(1,m'T) for 
m' > m. Let us choose I such that |L;| = 2"°/^|L|; in 
this case, because we consider H-dimensional systems, 
we have I = const • Then, the in¬ 

equality (35) reduces to (12). We thus prove Theorem 2. 


VI. ACKNOWLEDGEMENT 


We thank Yuji Shirai and Yoshihito Hotta for valu¬ 
able comments on the manuscript. TK acknowledges 
the support from JSPS grant no. 2611111. TM was 
supported by the JSPS Core-to-Core Program “Non¬ 
equilibrium dynamics of soft matter and information” 
and JSPS KAKENHI Grant No. 15K17718. K.S. was 
supported by MEXT grant no. 25103003. 

Note added. In updating the previous version 
(arXiv:1508.05797) to this version, related results in a 
different approach have appeared [115, 116]. 


Appendix A: Proof of Corollary 1 


The inequality (9) is shown by dividing as 

+ ||e-<“^^-e-*<’^||. (Al) 


The first term is bounded by (8). From Eq. (7), the sec¬ 


ond term is bounded from above by — Hp 

fin 


(")| 


T < 

< E:=„+i (;^(AT)-m!. For 
n -I- 1 < m < no, we have m! < n™“"“^(n -|- 1)!; hence, 
we obtain 




("O) 


T _ 


(ti)„ 


< Gln+lT 


in+2 


(A2) 


where we use l/(m -|- 1)^ < l/(n -|- 2)^ for rn > n -I- 1 
and XTuq < 1/16. By combining the inequality (Al) 
with (8) and (A2), we prove Corollary 1. ■ 


\\OL{t)-OL{t,i)\\ < \L\g{i,t), 
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Appendix B: Basic properties of multi-commutator 


For the analysis of few-body observables, we often 
treat the multi-commutators. Here, we show three lem¬ 
mas as useful technical tools in the analysis. For the 
convenience, we introduce two parameters {kA, Ja) to 
characterize basic properties of a generic operator A. 
We, in the following, refer to that an operator A is kA- 
local and J^i-extensive when 

A= ^ ax (/c^-local), 

\X\<kA 

\\ax\\<JA (J^-extensive), (Bl) 

X-.X3i 

for Vi G A. For example, the Hamiltonian given by (4) 
with (3) is fc-local and J-extensive. 

For the norm of multi-commutators, we can prove the 
following lemma: 

Lemma 3. Let {Ai}^^i be ici-local and -extensive, 
respectively, and Ol be an arbitrary operator supported 
in a region L. Then, the norm of the multi-commutator 

{A„, An-i,. ■. ,Ai, Ol\ '■= [A„, [An-i, , [Ai, O/,]] 

is bounded from above by 

n 

||{A„,A„_i,...,Ai,Ol}|| < X{{2JmKm)\\OL\l 

m—1 

(B2) 


where Km ■= \L\ + Y.i<Tn-i hi¬ 


proof of Lemma 3. We start from the following expres¬ 
sion for the multi-commutator [54]: 


{Ant ^n— 1 : ■ • ■ : Ai, Ol{ 


= E E E 

Xi\L X2\iL,Xi) X3\iL,XuX2) 


E 

X„\(L,XiX2...,X„-i) 


where we define Am '■= X)|x|<fc “ 

1,2, ...,n and X\{L, Xi, X 2 , ■ ■ ■, ^m) denotes the set 
of {A : A n (L U Ai U X 2 U • • • U A„) ^ 0}. We thus 
obtain 


[{Ati, An—It ...,Ai,0/j}|| 

E ■ ■ ■ E 11“^ 




Xi\L X„\{L,XiX2...,X„-i) 
For m = 1,2,3,..., n, we have 

E 

^ E 


(m) 


E 


Am) I 


(i/U-XiU-X2 U’•’U-Xtyi_ 1) Xjn'.XjTL^i 

< ^ Jrn 

(Z/U-Xx U-X2 U' ■ — 1) 

^Jm(\L\+ kij = KmJmt 

^ i<m-l ' 


where Jj„-extensiveness (Bl) implies the inequality 
Jm- We thus prove the lemma. 

As a related lemma, we also utilize the following 
lemma (See Ref [58] for the proof): 

Lemma 4. Let iL be a fc-local and J-extensive opera¬ 
tor. For an arbitrary fc^j-local operator A, the norm of 
the commutator [H, A] is bounded from above by 

||[Jf,A]|| <6JfcfcA||A||. (B3) 

The third lemma gives an upper bound of Ja for ar¬ 
bitrary multi-commutators: 

Lemma 5. Let {Ai}f_^ be fci-local with gi extensive¬ 
ness, respectively. Then, the extensiveness of the multi¬ 
commutator {A„, A„_i,..., Ai}, which we denote by 
A{a„,a„_i,...Ai}; is bounded from above by 

n 

A{21„,A„_i,...,Ai} ^ {‘^JmKm) (B4) 

m— 2 

with Km := J2i<m fc*- 


Proof of Lemma 5. 
induction. For n 

Sx Sx:ynX5^0[®x^ 
and hence obtain 


We 


prove 
2, we 




= EyE 


the lemma by 
have [A 2 ,Ai] 
( 2 ) 


the 


y Z^X:Xny#0L®x :“y J) 


(i)i 


J[A2Ad<E E 2||a^)||.||aW|| 

X3iy:ynX5.^0 

+ ^ ^ 2||a^)||.||«W|| 

y9ix:xny#0 

<2(A:i + fc2)JiJ2. (B5) 


We then assume that (B4) is true for n < uq — 1. 
For n = no, from (B5) and {A^^, A„o_i,..., Ai} = 
{A„o,{A„o_i,...,Ai}}, we have 


A{A„q,...,Ai} ^‘^JnoJ{Artg-i,...,Ai} [kng 
no 

= Jl \{{2JmKm)t 

m—2 


" 0-1 

E 


m—1 


where we use the fact that {A„(,_i,... A 2 , Ai} is at most 
fcm)-local. This completes the proof. ■ 

As a standard example, we apply the above lemmas 
to the case where |L| = /c, = fc and Ji = J for 

i = 1, 2,..., n. Then, Lemma 3 gives us 

||{A„,A„_i,...,Ai,Oi}|| < A^nlllOill, (B6) 


where we use the definition of A = 2k J. Moreover,from 
Lemma 5, we have 

J{A^,An-u■■■,A^} < A"-in! J. (B7) 

Finally, we denote that the FM expansion terms {fin} 
are fc^^^-local and jQ”^-extensive, respectively; from the 

inequality (B7) and Eq. (6), we have < (n -I- 1)A: 
and 


t(") 


< 


A"(n-H)!J 
(n -I- 1)2 


A"J 
n -I- 1 


nl. 


(B8) 
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Appendix C: Proof of Lemma 1 


Appendix D: Derivation of the inequalities (30) 
and (31) 


We here show the proof of Lemma 1, which gives the 
upper bound of {|lf^n||}5?Lo the FM-expansion (1). 

Proof of Lemma 1. For the proof, we start from the 
explicit form of the FM terms in Eq. (6). From the 
inequality (B6), We first obtain 


<{iL(t£r(n+l)): ■ ■ ■ j -f^(fo'(2))i ^(f-o-(l) ) } || 

+ \\{H{ta(n+l)), ■ ■ • ,^^0, i^(^CT(2))}|| 
<A"n! (Ikx(ia(i))|| + lkx(i<T( 2 ))||), 

|X|<fc 


We here evaluate the norms |ld)i_ji|| and which 

are defied in Eq. (28). It is enough to prove the inequal¬ 
ity for Ui and we can prove the other cases in 

the same way. 

We first calculate ||<i)i_„|j. From Eq. (21), it can be 
given by the expansion of 

= T[e~"-l'o (Dl) 

with respect to T, where we introduce the notation 

OO 

m—0 


which yields the upper bound of the norm Unnll as 




A”n! 

(n-f l)2r"+i 



E E (ll^x(Mi))|| + |kx(M2))ll) 


cr \X\<k 

2A^'n! .. / \ii 


2VbA’"n! 

(n+l)2’ 


where in the equality we use the definition of (5). This 
completes the proof of Lemma 1. ■ 


Each of the terms {rm(i)} is given by 



{Hi(tm), ■ ■ ■ Hi(ti), Vi(t)}, 

where Vi(t)} denotes the 

multi-commutator. By the use of (B6) with the defini¬ 
tion (20), we can obtain the inequality of 

||r™(f)|i<A- E lkx(t)||=:A-Vi(t). 

X:X 3 l.XnAi_i =0 

(D2) 

Note that Vi(t) < J from (4) and Y Vi{t)dt = TVi 
from the definition (26). Using the expression of (Dl), 
we have 


T"||4>i,„(r)|| <E E 

q>lmi-\ - \-mq—n—q'^^ ^ 


f'T nti r^q-1 


CH|rm,(u)l!triirm,(t2)|| • • •C’||r„,(t,)||dtidt2 


/O 


^ E E 

q>l miH- \-mq—n—q 


/ 7 ‘ 

/o Jo 


■ / Vi(ti)Vi(t2) ■ • • Vi{tq)dtidt2 ■ "dig 

Jo 

< E(AT)"2"-iEZE < Vir(2AT)"-^e^i/^ < ViTe^/(^'=)2-”+\ 


9>1 


(D3) 


where we use J2rm+-+m,=n-q = {Z-D ^^ the 
third inequality, e® — 1 < xe^ in the fourth inequality, 
and the condition T < 1/(4A) in the last inequality. We 
thus obtain the inequality (30). 


uY'^ = i for i = 1,2,..., TV; in the following, we con¬ 
sider no > 1. Remembering that vY'^ = — LI 

is given by the use of the FM expansion for H{f) and 
we first expand as 


Second, we calculate ||<I)i_„(T)||, which is given by the 


expansion of uY'^ = 


-A-n-o), 


r -i - '■dti 

T[e Jo 1 J 

with respect to T. For no = 0, we trivially obtain 

E OO 

n— 




n=no + l 






= 0 because of = Ho and 



E 

l,s—0 


(D4) 
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Then, we obtain a similar inequality to (D3): 


q>lmi-\-m 2 -\ - \-mq—n—q 

X iif:„ji-iir:„ji---r™ji, (d5) 


where we define ||f',„|| := 'Ei+s=m l|ri.s||- 

In Appendix D 1, we can prove the following inequal¬ 
ity for fiy. 


\\tiJ<{2Xy+^n‘oV^(^Y)- (D6) 


We thus obtain ||f^|| < Vi [2A(no-I-1)]< Vi(4Ano)'", 
where we use Uq > 1. The inequality (D5) then reduces 


to 

T"||4>("“)(T)|| <T"^ ^ (4Anor-«^ 

^ {SXnoTr [^i/(4Ano)]^ ^ 

~ ^ 9>1 ~ 

where we use T,m,+-+m,=n-q < 2””^ in the second 
inequality and the definition of uq = [l/(16Ar)J in the 
last inequality. This completes the proof of the inequal¬ 
ity (31). 

1. Proof of the inequality (D6) 

Remembering that H^°'^ and come from the 

FM expansion with respect to H{t) and respec¬ 

tively as in (22), we first expand 
as 

no no 

= (D7) 

q^O q^O 

where we denote the FM terms with resect to the Hamil- 

tonian Hi{t) by namely 

Using Eq. (6), the terms Rq is explicitly given by 


(9+1)2 


^(o')!(g ~ ^(o'))! 

_n ^ ^ 


il ,22 ,...ag + l —0,1 

{2l,22,...,2g + l}^{l,l,...,l} 


X ^qr^q+l (D8) 

where (t<„(,+i)), Mi^(tc,(,j)),..., (t^( 2 )). Mi,(fcr(i))} denotes the multi-commutator, Mo(t) := Vi(t) and 

Mi{t) := note that H{t) = Vi(t) + Hi{t) from the definition (20). 

We, in the following, calculate the expansion of e“*^i ^hyino)^ ih[ °h. 


OO OO 


{-tty 


s=0 /=0 q<l 


S\ 




(D9) 


91-1-92-1-I-93=*-9 


From the expansion (D9), we give the explicit form of 
in Eq- (D4) as 

h. = E E 

9<i 91-1-92-1- i-qs=l — q 

(DIO) 


In order to evaluate the norm of 
..., , i?g}||, we first derive the fol¬ 

lowing inequality: 


II {U;, M,,^, (Vi)i AF,, (f,),..., M,, (U)}|| 


< 


( 9 + l)A«i+«+i 9 i! 9 ! 
9i + 1 


Vl(ts-l-l) 


(Dll) 


for {ii,i 2 , ■ ■ ■ Jq+i} y {1,1,...,!}, where Vi(t) is de¬ 
fined in Eq. (D2) and we denote s such that ig+i = 0. 
For the proof, we define Ms '■= {Mi^(ts),... ,Mi,(<i)}, 
and obtain 


, ^iq + l (.^q+'^)^ ^iq (^ 9)1 ■ ■ • 1 (^ 1 )} 

~ {^91 ’ ^iq+1 (^9-l-l)’ ■ • ■ ’ ^iq+l (^s-l-l), ^s} 

= -{u;, M,,^, (Vi),... M., El(t.+i)}, (D12) 


where we use ig+i = 0 or Mi^_^^{ts+i) = Ei(ts-i-i) in the 
second equality. Because of the inequalities (B7) and 
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(B8), we have 

^91 ■ [(91 + l)A:]-local, -extensive, 

Mi{ti) : fc-local, J-extensive, 

Ms ■ (sfc)-local, [JA®“^s!]-extensive, 

and hence we obtain the inequality (Dll) by applying 
the inequality (B2) to (D12) with = X^Jnl/{n+l). 
From (D8) and (Dll), we have 


< Vi 


2'?+i(g-h l)A«i+9+i(?i!g! 

(<?i + !)(« + 1)^ 


(D13) 


where we use Vi{t)/Tdt = Vi and the coeffi¬ 
cient 29+1 comes from the summation with respect to 
■ ■ ■ ,* 9 + 1 } ^ in Eq. (D 8 ). In the 

same way, we can derive the inequality 


^ 29+iA9g! (g -H l)A9i + igi! {q + qi + 2)X‘i^ + ^q 2 l 

“ ^ (g-l-1)2 gi -I-1 92 + 1 

(9 + 9i + 92 + • • • + 9s-i + s)A9=+1(7s! 

9s + 1 

<Vi29+iA9+"9!A9i+9^+ '+9»qJq,_i! • • • gi! 

(9 + 9i + 2) • ■ ■ (9 + 9i + 92 + • ■ ■ + 9s-i + s) 
<Vi29+iA'+"n[,(Z -h 2) • • • (/ -h s), (D14) 


where in the last inequality we use the conditions 9 < 

no, Qj < no (j = 1 , 2 ,..., s) and 91 + 92 +-h 9s = 

I — q, which yield the inequalities 9!9s!9s_i! • • • 91! < 

- +91 ^ -< 

l+j for j = 1,2,... ,s. 

From Eq. (DIO) and the inequality (D 14 ), we finally 
obtain the inequality (D6) as 

WhsW < ^ 29+i2'+^-9-i(Z + 2) • • • (Z + s) 

q<l 

= ^!^^n[,2'+«(/ + 1)(; + 2 )---il + s) 

Si 

={ 2 xy+%v,(^Y^, 

where we use E9i+92+-+9.=i-9 = < 

2l+s-q-l 

in the first inequality and J 2 q<i — ^ + 1 in 
the second equality, respectively. 


Appendix E: Proof of Theorem 3 


The theorem gives the upper bound of the norm 
Then, we start from the inequal¬ 
ity” 


l|n 


(n) 


U^U 


(n) 


>E+AE'^F ^^<E 


^<E\ 




-rH 


("o) , 


tH, 


{") 

F e 


-rH, 


("o), 




(^) , 


(El) 


where we use W < and 

||grffp < e'^^. We therefore evaluate the 

three norms: ||e'^'^F''^ ||, ||e'^-^F ' || and 

IlgFffp {7pe“'^^F II. We first show the sketch and give 
the technical details in Appendixes E 1 and E2. 


First, to evaluate the norm of 
use 


r,rH, 


("o) 


-tH- 


(") , 


we 


„-rH, 


(") 


= e 


-tH, 


("O) 


r 


r 

eJo 


F F dx 


which yields 


|e"<“>e-"^. 


(») , 


< 


/ pT TLq \ 

exp(/ £ T™||e"<“’fli„e-"<“’||cix). (E2) 


By using Lemma 1 on the norm of we obtain the 
upper bound of 

||e"<’e-"<“'|| < gFVoT"+B+„+i/2^ (E3) 


with the derivation in Appendix E 1, where W 71+1 is de¬ 
fined in the theorem: Wn = 2(4A/3)"'n!/(n -|- 1 )^. For 
||grff^ ^ 11 ^ ggg g^gQ Q] 3 |;a,in the same inequal¬ 

ity. 

We then consider the norm of 


„rH, 


("o). 


C/pe 


-rH 


("O) , 


< 


„rH 


("O) 


UfU, 




where we use [C/p”°\7Lp”“^] = S, ||C/p”“^|| = 1 and 
the inequality ||Ai3|| < ||A|| • ||i3|| for arbitrary oper¬ 
ators A and B. We here adopt the decompositions 
(18) and (19), and introduce the notations of Si = 

||e"<”VF.i(T)C/^”“^^(T)e-<‘”|| for i = 0,1,...,iV. 
Because of [/F,i-i('r) = UF,i{T)Ui from Eq. (21), we 
have 

si-i < ||e"<°VF,i(T)C/^7>^(T)e-"<“’|| 

+ ||e"<‘’VF,i(r)(l - c/it7f“^^)[/^7^^(r)e-"<”’|| 
< Si + 44Vi JT22""«/^ (E4) 


where the second inequality comes from lemma 2 which 
ensures C/iC/f"°^^ ~ 1; we give the derivation in Ap¬ 
pendix E2. From the definition of Si, we have sq = 
||gFffp °’jjpf 7 +o)tg-Ti/(, o’li gggj ^ Hence, from 

Eq. (27), the inequality (E4) reduces to, 

||e^4’'”VFe"^-f^F"“’|| < So < 1 + 22TVbAT2-"»/2^ 

(E5) 


where we use the assumption T <t and J = A/(2fc) < 
A/2. 

By applying the inequalities (E3) and (E5) to the 
inequality (El), we prove Theorem 3. I 
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1. Derivation of the inequality (E3) 


in the similar way to the derivation of (A2). 


We start from the inequality (E2) and evaluate the 

norm of ||' || • For the purpose, we utilize 
the following technical lemma: 

Lemma 6. Let A be an arbitrary (nAk)-local opera¬ 
tor such that A = Then, the norm of 

o>|| jg bounded from above by: 




(r^o) 


Ae 


-xH, 


("O) , 


< 


2(l-2r;AnoT) 


Pll, (E6) 


where ry := 1/(1 — 2Aa;) and A := 6k^J. The proof is 
given below. 


By combining the inequalities (E2) and (E7), we ob¬ 
tain 'll < same in¬ 
equality can be given for ||e'^'^F ° ||. We thus 

prove the inequality in (E3). 


Proof of Lemma 6. We here calculate an upper bound 


of lle^-^F 


("O), 


Ae """F lb where A is an arbitrary (n^A:)- 
local operator due to the condition in Lemma 6. We 
have 


From Eq. (6), the operator Llm is given by the sum 
of (to -|- l)th multi-commutators, and hence it contains 
at most (to -I- l)fc-body couplings, namely ua = m+ 1. 
For X < T = 1/(8A) and no = [1/(16AT)J, we have 
f] < 4/3 and 1/[2(1 — 2r]XnoT)] < 3/5. Therefore, using 
Lemma 6 with Lemma 1, we obtain 


no 

E 


rpn 




(r^0) 


Llr 


^—xH, 




<! y; 

5 

7n—n-\-l 


if)-.. 


^VoiUT/3r+\ , EoT"+iW„+i 

- (n + 2)2 - 2 


(E7) 


E 

mi+7712 + ...+ 7713 

{mj<no} 


1 , ^nis — i 7 • ■ • 7^mi7A|-||. 


Under the assumption of toi -I- TO 2 -I- • • • -I- ms = M, we 
use the inequality (B3) and obtain 


||{f^ms;^ms-i7 ■ • ■ 7 7^/11 

s 

< JJ[6JQ™''^A:^(TOj -I- l)]n^(nA + mi + l){nA + mi + m 2 + 2) ■ ■ ■ {ua + mi+m 2 + ■ ■ ■ + mg-i -I- s — 1)||A|| 
i=i 


<(6Jfc2)«A™i+™=+-+™= toi!to 2 ! • • • ms^UAinA -I- toi -|- 1) • • • {ua + mi + m 2 + ■ ■ • + mg-i -I- s — 1)||A|| 
<W(Ano)^(nA + M)(nA + M + 1) • • • (n^ + M + s - 1)||A||, (E8) 


where we use the fact that Llm is [(to-I- l)A:]-local and jQ™^-extensive with < X^Jm\/{m+ 1) as in (B7), and 
to! < n™ for to < no- We thus obtain the inequality (E6) as follows: 




(^ 0 ), 


00 00 


< ll^ll E E eiA eE2lE(„, + M){n^ + M + 1)... („.4 + M + . - 1) 


S! 

M=0 S=0 

(1 - 2Ax)-"'^ 


E 

M=0 


2XnoT 

l-2Xx 


M 


2(1 - 2rjXnoT) 


II+I7 


where the first inequality comes from (E8) and 

equality comes from the Taylor expansion of (1 — 
2Aa:)“^"'^+'^!* and we use the definition of 77 := 1/(1 — 
2Aa:) in the second equality. This completes the proof 
of the lemma. I 


2. Derivation of the inequality (E4) 

We here derive the upper bound of Ai := 
||e"<‘’’c/F.*(T)(C/,C/f“)^ - i)C/^^^,(T)e-"<“'|| as A, < 
44Vi JT^2“"“/^. For this purpose, we first expand 
respect to T and obtain 

00 00 

m,m'—0 n—0 
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where we use the notations in (28) and define := 

Because is equal to 

UiUj up to the order of T"o, we have 

^ (E9) 

n>no + l 

and hence we obtain 

A, < ^ (ElO) 

n>no + l 


and evaluate the norm of ’^i^n+q- From the defini¬ 
tion (21) of U-p^i{T), we obtain 




n+q 





Because Hi{ti) is fc-local and J-extensive, the in¬ 
equality (B3) gives 


where A,,„ := VF.*(r)^,,„4_,(T)e-"^ 

We prove the following inequality below: 



Kn < 

5 


(Ell) 


< \'^n{n+ l)---{n + q - l)|j^'j,„||, 

where A = and use the fact that is at most 
(nfc)-local. We thus obtain 


for T = 1/(8A). Moreover, by following the derivation 
of the inequalities (30) and (31) with T < 1/(8A) < 
1/(8A), we have 

T™||$i,m(r)|| < v,Tei/(2fc)2-2m+2^ 


m^^,n+q\\ < (Ar)^ ^(^+^)--> + g 

(E14) 

By substituting A with 'i'i^n+q in the inequality (E6) of 
Lemma 6, we have 


which yields 


—n 

< 15VfT^ 2-2™-'"' < 30ViJr22-"^ (e12) 




,n+qt 


„-rH: 


("o) , 


- 2(1 - 2r]XnoT) 


(E15) 


with rj = 1/(1 —2Ar). Thus, from the inequalities (E14) 
and (E15), we obtain 


where we use 8e®/® < 15 and Vj < J. 

By combining the inequalities (Ell) and (E12) with 
(ElO), we obtain 

A* Y < 18ViJT2 Y 2■"''^ 

n>no + l n>no + l 

which finally reduces to the inequality of A^ < 

44ViJr22-"«/2, 

Derivation of the inequality (Ell). For the evaluation 
of Ai n-! we first expand as 


C/F.*(T)4',,„C/if ,(T) = ^ T^^,^n+q (E13) 

9=0 


A, = 




(^o) 


?7f.*(T) 4',(r)e 


-tH, 


("o) , 


{\TYrf-+‘i n{n+l)---{n + q-l) 
- ^ 2(1 - 277AnoT) q\ 


2{l-2ijXnoT) \l-r]XT 


(E16) 


where we use the fact that ^i,n +9 is [k{n + g)]-local 
in the second inequality and the last equality comes 
from the Taylor expansion of (1 — rjXT)~'^ with respect 
to -qXT. We here set r = 1/(8A), no = [1/(16AT)J 
and T < 1/(8A). Hence, we have 77 = 4/3, Y^xt — 

= I < 2 I /2 and < 3/5, from which 

the inequality (E16) reduces to (Ell). 
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